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Abstract 

We study the question of existence and uniqueness for the finite 
temperature Kohn-Sham equations. For finite volumes, a unique 
soluion is shown to exists if the effective potential satisfies a set of 
general conditions and the coupling constant is smaller than a certain 
value. For periodic background potentials, this value is proven to be 
volume independent. In this case, the finite volume solutions are 
shown to converge as the thermodynamic limit is considered. The 
local density approximation is shown to satisfy the general conditions 
mentioned above. 

Key words: density functional theory, Kohn-Sham equations, ex- 
istence and uniqueness, thermodynamic limit, periodic potentials. 
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1 Introduction 



For a self-interacting TV-body quantum system with a background potential 
m, the finite temperature Kohn-Sham equations consist of: 



The exchange-correlation potential v xc [n] may have a complicated func- 
tional dependence on the particle density. However, it acts on the Kohn- 
Sham orbitals as a regular potential, by multiplying the orbitals with a 
function. yQ 

These equations determine the density of particles at the thermody- 
namic equilibrium of the system. The zero temperature formalism proposed 
by Hohenberg and Kohn 2 , now known as the Density Functional The- 
ory, can be extended to finite temperatures if, instead of the ground state 
energy, one considers the grand canonical potential. |3] In this case, one 
can show that the grand canonical potential is a functional of the particle 
density. This functional achieves its minimum for the equilibrium particle 
density. Eqs. Q |2J represents the Euler-Lagrange equations associated 
with this functional. |S] In contrast to the zero temperature case, the 
finite temperature Kohn-Sham equations involve an infinite set of orbitals. 
In practical applications however, one only has to consider the orbitals with 
energies up to /j + ksT. Above this limit, the contribution to the density of 
particles Eq. (J2J) becomes negligible. The finite temperature formalism has 
a major advantage over the zero temperature formalism in that it avoids 
the problem associated with the degeneracy of the last occupied energy 
level. 

While most of the work on these equations has been concerned with 
rigorous derivations or finding better and better approximation of the 
exchange-correlation potential, little it is known about their solutions. |SJEI 
Computational physicists often assume that the Kohn-Sham equations have 
a unique solution. The same assumption has been made for long time for 
models like Hartree or Hartree-Fock. Despite many efforts however, the 
uniqueness is still an open problem for these models [Hj and only little ad- 
vance nas been made in this direction. Moreover, the symmetry 
breaking within these models |12l IT5| shows that the uniqueness can be a 
very delicate problem. 

The first goal of this paper is to search for those conditions on the 
effective potential that guarantee the existence and uniqueness of a solution 




(1) 



where the value of the chemical potential is found from: 
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for the Kohn-Sham equations on a finite volume. The conditions will be 
formulated at an abstract level, without making references to any explicit 
expression or approximation of the exchange correlation potential. We will 
show in the last section that these abstract conditions are general enough 
to include the local density approximation. 

The Kohn-Sham equations are derived from the Kohn-Hohcnberg 2 
functional. This functional is based on the assumption that the density 
of particles is u-representable. Despite many efforts, the u-representability 
problem has not been yet rigorously solved for infinite volume. A functional 
that does not require w-representability has been proposed by Lieb. m] Un- 
fortunately, this new functional does not automatically lead to the Kohn- 
Sham equations. It was shown however that, for finite volumes, the Kohn- 
Hohcnberg and Lieb functionals coincide. £Q The relation between the two 
functionals is not yet completely understood for infinite volume. The ther- 
modynamic limit of the Kohn-Sham equations can therefore be considered 
as a fundamental issue in density functional theory. We believe that the 
unsolved part of the w-representability problem can be avoided by study- 
ing the thermodynamic limit of the finite volume Kohn-Sham equations 
instead of considering the infinite volume. 

The second goal of this paper is to find the general conditions on the 
effective potential which guarantee that the finite volume solutions of the 
Kohn-Sham equations have a well defined thermodynamic limit. Our solu- 
tion to this problem applies only for periodic background potentials. 

Within the local density approximation, partial results on the Kohn- 
Sham equations have been reported in Ref. JHj- The methods devel- 
oped there could not be used for realistic exchange-correlations potentials. 
The main obstacle was the low density behavior of realistic exchange- 
correlations potentials. In the last section of this paper we will show how 
this problem can be solved. We will show that the local density approxima- 
tion of the exchange-correlation potential satisfies the abstract conditions 
mentioned above and a unique solution exists. 



2 The fixed point approach 

We assume in the following that the background potential comes from the 
interaction of the particles with background charges of opposite sign, which 
are considered fixed and given. In this case, the Kohn-Sham equations take 
the following form: 

J (-±A + (n-n )*v + v xc [n])r) m =s m r) m 
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We also consider a system which is charge neutral: 



(n (x) - n (x)) dx = 0. (4) 

vol 

The chemical potential must be adjusted such that the above condition is 
satisfied. The neutrality condition will play an essential role when long 
range interactions are considered and it will help us to improve our pre- 
vious estimates on the Hartree potential. We borrowed the idea from 
the homogeneous electron gas where it is known that the neutrality con- 
dition is essential when the thermodynamic limit is considered. Also, 
the neutrality or partial neutrality condition play an essential role in the 
Hartree-Fock model when applied to Coulomb systems. ^3 El The neu- 
trality condition will be further discussed at the end of Section 3. 

Let us denote the Kohn-Sham Hamiltonian and the effective potential 

by: 

H n = -±A + AV>] (5) 
V [n] = (n-n Q )*v + v xc [n] , 

and the Fermi-Dirac distribution by: 

$FD(t)=(l + eP t )~\ (6) 

For finite volumes, we now formulate the Kohn-Sham equations as a fixed 
point problem. In this case, A is the Laplace operator over the volume vol 
constrained by various boundary conditions. 



Theorem 1 (The fixed point approach) Let S N C L 1 (vol) be defined 
as: 

S N = {n e L 1 (vol) , \\n\\ L1{vol) = N\ . (7) 

Suppose that for n E S N and a > the following condition is satisfied: 

V>](-iA + a) -1 | < 7a , (8) 

where j a may depend on N. Then the map: 

T : S N ^S N (9) 
S N 3 n -> T[n] (x) = <S> FD (H n - /i„) (x, x) 

is well defined. Here, fj, n represents the unique solution of the equation: 

N = Tr$ FD (H n -fi n ). (10) 

Moreover, the fixed points of the map T generates all the solutions of the 
Kohn-Sham equations. 
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Let us prove first the following result which will be used many times in 
the following. 



Lemma 2 Let Hq be a self-adjoint, bounded from below Hamiltonian and 
suppose that exp (—Hq) is of trace class. Let V be a self-adjoint potential 
such that, for a > |inf a (-ffo)l : 



V(H + a) 1 s: 7o . 
Then H = Hq + XV is self-adjoint and: 

Trcj>^ [(I + A7JH0 + A07J 
sc: Tr (l + e^ H -^y j 



(11) 



(12) 



^ e 



jf}{fi+\a'Y a 



Tr e 



-jf}{l-X la )H 



where 4> J ■ , j = 1, 2,..., arc monotone decreasing, convex functions on 



[eo, 00), eo = infer {H), such that: 



1 + e 



/3(*-M) 



All functions <pj „ (f ) can be chosen such that lim^^ 
(*) > 0. 



(*) 



(13) 



1 and 



Proof. We will use classical techniques from Ref. ^Sl ■ We start with the 
first inequality. Let {iJmj £ m} m be the eigenvectors and the corresponding 
eigenvalues of Hq. rf m ^V (H) and we can write: 

Tr (l + e^ H -^Y 3 (14) 

> Tr^(H)^^2 m ^((r,° m ,Hr,° m )) 

= E m h» + a ) (i + W (Ho + ay 1 ) rt) - «) 

> H m ^(( 1 + A 7a)£m + Aa7a) 
= Tr^((l + A 7o )ifo + Aa7j. 

For the second inequality, let {j? m) E m } be the eigenvectors and the cor- 
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responding eigenvalues of H. Then: 

= e J/3(a+M)y^ e -i/3(l->7j(e m +a)( J ) m ,(/+Ay(Ho+a)- 1 ) _1 7 ?m ) 

J> e i/3(a+M) \^ g-J/3(l-A7 )(£m+a)/(l-A7 <> ) 
— 'm 
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Proof of Theorem 1. We need to show that Eq. (|10fl has a unique 
solution for all n € S . From the previous Lemma, 

-j—Tr $fd (H n — fj,) (16) 
a/i 

> 0e-^Trcj> 2 ^ (-I(1 + A 7 J A + Aa 7o ) . 

The last term is strictly positive when we choose 2 u > 0- Also, from the 
previous Lemma, it follows that the above derivative is finite. Then the 
right hand side of Eq. 1(10(1 is a strictly monotone, continuous function of 
fi. As we already mentioned, </>,•„ can be chosen such that <fij „ — > 1 as 
goes to +00. Then one can see from 

Tr^[-±(l + \ la )A + \a la \ (17) 
< Tr$ F1J (ir n -/i) 

^ e 0{H+\aj a )j, re f3/2(l-\ la )A 

that, as is varied from — 00 to +00, the right hand side of Eq. 110(1 
varies from to +00. This, combined with the strict monotonicity and 
continuity, allows us to conclude that Eq. 1 (10(1 has a unique solution. 
Using the eigenvectors and the eigenvalues of H n , 

(-5A + (n - n ) * v + v xc [n]) r\ m = £ m r) m , (18) 

the fixed point equation for T reduces to 

n (a?) = (l + eP^-"J) ~ l \ Vm (x)\ 2 . (19) 

m 

Eqs. ((18(1 and ((19(1 represent exactly the Kohn-Sham equations. ■ 

The fixed points of the map T generates all the solutions of the Kohn- 
Sham equations because S N is the largest set where these solutions can be 
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found. Another important observation is that there exists an upper and 
lower limit on the chemical potential, limits which may depend in general 
on the number of particles. This can be seen from: 



Trt t>i# n + A 7a) A + Xa la ] (20) 

TV < e / 3 (M„+Aa7„)y re /3/2(l-A 7 JA^ 



For periodic background potentials however, we will show that this limits 
are independent of the number of particles. 

3 The Kohn-Sham equations with periodic 
background potentials 

In the following we consider a background charge distribution which is 
periodic with respect to a lattice T, 

r = jxe R 3 , x = Y^i l 5 ll n l £ zj , (21) 

i.e. no (x + Rj = n a (x) almost everywhere when fie T. 5 i7 i = 1,3, 

represent three linearly independent vectors in R 3 . Let us also consider a 
finite crystal confined in the volume: 

V = |f e R 3 , x = sC a* < if | , (22) 

where if is a positive integer. We denote the crystal's unit cell by: 

cell = \ x e i? 3 , x = \ ' a'Si, sC a 1 sC 1 } . (23) 



i=i J 



Thus, the crystal is formed from K 3 unit cells. We impose periodic bound- 
ary conditions and we also allow the particles on opposite faces of the 
crystal to interact to each other. The resulting problem is that of parti- 
cles trapped on a torus T obtained by connecting the opposite faces of the 
crystal. Any point from R 3 can be viewed as a point of the torus. The 
kinetic term of the Kohn-Sham Hamiltonian is given by — |A, where A 
represents the Laplace operator over the torus T. We will assume that 
the particles interact via a two-body potential which depends only on the 
distance between particles: 

v(x,y) = v{\x,y\) , x,y£T, (24) 
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where |-,-| denotes the distance on the torus. In this case, the potential 
generated by the background charge is T-periodic. Indeed, for R £ T: 

v (^x + R,y^jn (y)dy = J v ( x,y - R^j n (y) dy (25) 

v(\x,y\)n (y + Rj dy 
v(\x,y\)n {y)dy, 



where we also used that the measure dy is invariant at translations. We 
will assume that the exchange-correlation potential is T-periodic when the 
density of particles is T-periodic. In this case, the set of T-periodic density 
of particles: 



S^ er = {ne S N ,n(x + R\ =n{x) a.e., Re r} 



(26) 



is invariant for the map T. In this paper, we will search for the fixed points 
of the map T only in this invariant set. Thus, from now on, we will restrict 
T to Sp er . For finite volume, the above system include also the case of 
nonperiodic background potentials. To include such systems, the unit cell 
is taken equal to the entire. Let us denote by iVo = N/K 3 the number of 
particles per unit cell. The thermodynamic limit is defined by fixing No 
and letting the number of unit cells to go to infinity. 
Let us consider the following unitary transformation: 



L 2 [cell] 

qeAjr 



U : L 2 (T) 

q 

L 2 (T)3f^ (Uf), 



(27) 



raGAj. 



where 



2?r 
~K 



q and A K = {0, 1, K - l} 3 . For n e S. 



N 

per > 



Sham Hamiltonian is T-periodic and consequently: [T^] 



UH n U-'= {-\* K + \V[n]) = H£\ 

qeAjr qeA K 



the Kohn- 



(28) 



where is the Laplace operator over the unit cell with the following 
boundary conditions: 



/ (f + Sj)= e w3 f (x) and /' (x + 6,) = e 103 f (x) 



(29) 
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for x and x + 8j on the faces of the unit cell. The symbol /' stands for 
the derivative of / along 5j. V [n] in Eq. (|28[) is just the restriction of 
the effective potential to the unit cell. Because the kinetic energy and the 
effective potential depends on the volume, we will write the Kohn-Sham 
Hamiltonian as: 

=-\& + V {K) [n\. (30) 

We will write Tk to indicate that the map defined in the previous section 
depends on K. Also, because N = K 3 Nq, it will be more convenient to use 
the notation S^ er instead of Sp er . 

Let us discuss now the neutrality condition. This condition seems ar- 
tificial for interactions other than the Coulomb force. However, when the 
volume is transformed to a torus and the particles interact via a potential 
which depends only on the distance between the particles, one can imme- 
diately see that adding a uniform background charge has the effect of a 
constant added to the Kohn-Sham Hamiltonian or to the full many-body 
Hamiltonian. This does not affect the solutions of the Kohn-Sham equa- 
tions or the physiscs of the problem. Thus, we can allways add a uniform 
background charge such that the neutrality condition Eq. (£IJ is satisfied. 
The neutrality condition can be regarded as a mathematical artifact. 



4 The result 

Our abstract conditions for existence, uniqueness and thermodynamic limit 
for the Kohn-Sham equations consist of the following. 

(CI) For any 9 G [0, 2tt) 3 , n G S* r and a > 0, 



(31) 



where it is assumed that 7 Q depends only on a (when Nq is kept 
fixed). 

(C2) There exists a closed set B C L 1 {cell) such that T K [Sf er \ C B K , 
where 

B K = {n G S£ r , n\ cell G B) , (32) 

and, for n.1,2 € Bk, there exists a constant L, independent of K, such 
that: 



V {K) [ni] - V {K) [n 2 



^L\\ni-n 2 \\ L i icell) . (33) 
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(C3) Forne 



per ' 

y{K+i) 



as if — > oo, 



(34) 



where n in yC*^ 1 ) [71] represents the unique extension of n in S^ 1 . 
Our main result is given below. 



Theorem 3 Suppose that (Cl)-(CS) are satisfied. Then: 

i) The maps Tk are well defined. 

ii) Tk have a unique fixed point provided the coupling constant is smaller 
than a certain value which is independent of K. 

in) The thermodynamic limit of the fixed points is well defined. 

Because the proof of the first two points and the proof of the third point 
of the above Theorem are based on different techniques, we present them 
in two separate subsections. 

4.1 Existence and uniqueness for finite volume 

Proof of Theorem 3. i) Using the unitary transformation Eq. (J57J, for 
/6l 2 (T): 




< la E ll/qlli 2 (ceH) =7oll/lli»(T)- 

Thus, the conditions of Theorem 1 are uniformly satisfied and in conse- 
quence, all the maps Tk are well defined. From this uniform estimate one 
can easily conclude that the spectrum of H^'^ is uniformly bounded from 
below, i.e. there exists an eo, independent of K, such that a 

[eo, 00) and, consequently, a (H { n K) ^j C [e , 00) for all K. We can also show 

that the upper and lower limits of the chemical potential do not depend on 
the volume. Indeed, using again the unitary transformation Eq. I|27|) . the 
inequalities Eq. I|2U|I can be transformed into: 

inf Tr^ (-§ (1 + A 7a ) A 9 + Aa 7a ) (36) 

< No < e^ +Xa ^ sup Tre^ 1 -^)^, 
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where the infimum and supremum goes over all 6 e [0,2tt) . Using the 
explicit expressions for the eigenvalues of A^ and the fact that we can chose 
<p 1 such that lim^oo <j) 1 — 1, we can conclude from above that there 
exist \x m and /i M , independent of K, such that /x„ G \p m , M M ]. For A 1, 
all these parameters, eo, Mm an d fMi can be considered A independent. ■ 
We prepare now for the proof of the second point of Theorem 3. Let us 
prove an estimate which will be used many times in the following. Along 
this paper, 1 1 - 1 1 ^ will denote the trace norm. 

Proposition 4 Let F be an analytic function in a vicinity of [eo, oo) and 
V\ and Vi two self-adjoint potentials over the unit cell such that: 



vi, 2 H&e+oy 1 



Then, 



\\F (-±A S + XVi) F (-|Aj + AV2) I, 
< ^7f \\Vi-V 2 



(37) 



(38) 



2Nii(ceH) ' 

where j F depends only on the function F. 

Proof. Using the notation AU = V\ — V<i and g a:Z (x) = (x + a) / (x — z), 
we can write after simple manipulations: 



A 



F(-lA d + \V 1 )-F(-lA d + \V 2 ) = 
x J dzF(z)g a , z (— + AVi) (l + A (-§A ? + a)" 1 Ki) 

x(-5 A e - + a ) _lA ^(-^ A 5 + «) _1 

x (l + AU 2 (-iAg + a)" 1 ) ~* g a , z (-|A ? + AU 2 ) , 



(39) 



where the integral is on a curve that surrounds [eo, 00) and belongs to the 
analyticity domain of F. In consequence: 



(40) 



Let us consider: 



A 



(1 - A 7 J' 



■ / \ dz \ \ F ( Z )\ SU P \9a,z [X 



AV 



a;el£o ,00) 

IV 



H A * + < 
A=\AV\ 1/2 i-^ + a)- 1 



(41) 
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where the square root is defined through the polar decomposition: AV 
S |AV|. We can immediately see that A is a Hilbert-Schmidt operator: 



\A*A\ 



dx 



cell 



AV(x)(-^A s + a) 2 {x,x) < k a ||AV|| il(ceM) . (42) 



We used the fact that (— |Ag + a) 2 {x, x) can be computed explicitly 
and: Q3] 

(~^A s + a y 2 (x.x) ^ k a , (43) 
with k a independent of 9 £ [0,27r) 3 . Then we can continue: 



- WA^ SA^ ^ \\A f \\ HS \\SA\\ HS 

< \\M 2 hs = \\^a\\, 



(44) 



and this, together with Eqs. I|40(l and l|42|l . proves the affirmation. We can 
also identify j F : 



If 



dz\ \F(z)\ sup \g a ,z{x)[ 

x£[e ,oo) 



(i - a 7q ) 

Proof of Theorem 3. ii). For m,2 £ Bk one has: 
||T[m] - r[n 2 ]|| L i (r) 



(45) 



(46) 



<f>FD ( Hg > - M ni 



(•^0 _ u 

n 2 ™ 2 



Using the monotonicity of the Fermi-Dirac statistics with respect to the 
chemical potential we can write: 



(47) 



Mm 



At this point we use the fact that: 



(48) 
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so we can conclude: 



\T[ni] -T[n 2 ]\\ L i [T) 



(49) 



s$ 2 



\*FD (HW-(l n2 



'FD 



JJ ni 



- ®FD 
Vn 2 



= 2 E I* 

< 2AK 3 7$ |yW[m]-yW[n 2 ] 
sC 2\j$L || m - ri2|li,i(T) • 



H-(K,q) 

ri2 



A simple analysis of the expression l|45|) shows that 7$ is maximum when 
/i„ 2 = 7$ can be also chosen A independent for A ^ 1. Thus we 

proved that, for A smaller than a certain constant, independent of K, the 
maps Tk are contractions on the closed, invariant sets Bk- This implies 
that they have a unique fixed point in Bk- Because Tk [Sp er ] C Bk, it 
follows that Tk have a unique fixed point over the entire S^ er M 
We end this section with estimates on the chemical potential. 

Proposition 5 For [i x 2 S [/j m ,/i M ] and n S Sp er , there exists C and C 
strictly positive constants, independent of K such that: 



C K - M2I < K- 3 \F K (n, Ml ) - F K (n, M 2 )K C" K - M2I 



where: 



F K (n, M ) =Tr$ FD 
Proof. The affirmation follows from: 



I-Pk (n, Mi) - (n,/i 2 )| 



, dF K (n >f i) 



/'1 



<3/i 



and from estimates on 

dF K {n, n) 



1 + e 



/3(^i K) -M) 



These estimates follows from: 



p e -0»MTr (l + e 



(50) 
(51) 

(52) 

(53) 
(54) 
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which can be reduced to: 

Pe-to" inf Tr^ m (-1 (1 + A 7 J A ? + Aa 7a ) (55) 
s; K -* dF *( n 'ri <; pe-P^JPQtM+iWa) sup Tr e^^^^s, 

by using the unitary transformation i)27|) and Lemma 2. The infimum and 
supremum go over all 9 G [0,2tt) 3 . We can conclude that: 



C < K -* dF *} n >ri < C > 



(56) 



where C is strictly positive if (f> 2 > 0, C < oo and both constants do not 
depend on K. For A ^ 1, these constants can be chosen independently of 
A.I 

Theorem 6 For any n G Bj<, the sequence s/iT° m [n] r converges to a 
unique limit. 



Proof. For 711,2 € -Br-, we have successively: 

C# 3 |/V - /in,] 

|F K (7ii,/J ni ) - F K (m,M na )| 

= |iv-i^(m,/i n2 )| 

= l-PjST (n 2 ,Mn 2 ) ( n l'Mn a )| 



(57) 



FZ3 



AT) 



y ni 



) _ 



Aif 3 7 $ L ||tii -712 11^ 



The affirmation follows from the fact that {T om [n]} converges to the 
same limit for any n £ S, 



K 

per ' 



4.2 The thermodynamic limit 

For a given K, let us denote by uk and /i^ the fixed point and the corre- 
sponding chemical potential of the map Tk ■ We prepare now to prove the 
last point of Theorem 3. We will use the following result from Ref. |11| . 

Proposition 7 Let fi 6 IPthi^m] an d ^ be a self-adjoint potential such 
that: 

"v(-±A d + a y l ^ 7a , (58) 
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for all 6 e [0,271-)"*. Then: 

\\9 FD (-%A 3 + V-fi) -$f/3 (-sA r + V-fi)\\ 1 < ct. 
where ci. and e depend only on j a . 



(59) 



The thermodynamic limit of the Kohn-Sam equations will follow from 
the following result. 



Lemma 8 With the above notations, 

\\n K+ i - n#|| il(ceM) -> as X -> oo. 



(60) 



Proof. There is a unique extension of tlk in which will be denoted 

by the same symbol n^- An important observation is that nx € B^+i- 
Let us prove first that 



( H (K+1) 



UK 



(x,x) 



(61) 



goes to zero as K goes to infinity. Indeed, Eq. (|61|l can be evaluated as it 
follows: 



w E * M (fl^ +1|q) -/*jt)(Sf,50- 

qgA A -+l 



(62) 



qeAif 



L 1 (ce/Z) 

qG3A A -+i 



h7 E |M^« +1,q) -^)-M^ ,q) 



qeA* 



where c?Ax+i = A^ + i\A#-. The first two terms above go to zero as K goes 
to infinity. For the last term we use the following: 

\\Q FD (-|A ? + V - n) - t> FD (-iA r (63) 
s? ||$ FD (-§Ag + V - m) - *fd (-|A g + V - /x) ||j 
+ ||3 FD (-±A e - + V' - Mif) - (-|A r + V - m) || x 

< A7*||V-VV(od0 +rf - 
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Replacing and 9 by 27rq/ {K + l) and 2nq/K, V and V by [n*-] 
and V {K) [uk] and /i by \i K , Eq. (J^TJ follows immediately from condition 
((73). 

Let us denote by fi nK be the chemical potential corresponding to uk 
when riK is extended in Bk+i, i.e.: 

(X + l) 3 No = Tr $ FD - 



It follows that: 

Indeed, from Proposition 5, 



as K 



■ oo. 



(64) 
(65) 
(66) 



Tr<j> 



(^ +1) -MA-)-rr$ 



®FD 



\ n 



K+l) 



Then Eq. (JBSJl follows from Eq. igTJ. 
Finally, 



OO 

Efei M-T°tr x) [»*] 



(67) 



L!(ceH) 

< (1 - 2Ai 7tI> )" 1 ||T A '+i K] - ^if|lz,i( ceH ) 



and: 



\\T K +i [nic] ~ n K \\ L1{cell) 



(68) 



n k 



(x,x) 



II K 



t 1 (ce(i) 



L 1 (ceH) 



L 1 (ceZ!) 



We already proved that the last term goes to zero as X — > oo. Using the 
monotonicity of the Fermi-Dirac, the second last term of Eq. (|68|l is equal 
to 



(K + l) 3 (n K ,/j, nK ) - ("if,^)) «S C" |^ njl 



Mid 



(69) 
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where we used the same notations as in Proposition 5.1 

Within our conditions, we cannot prove that {uk}k is a Cauchy se- 
quence and in consequence the problem of the thermodynamic limit is not 
yet solved. However, using the weaker result of the above Lemma, we can 
prove that the density of particle converges in a distributional sense as the 
thermodynamic limit is considered. This will end the proof of Theorem 3. 

Theorem 9 If tik is viewed as a linear functional over L°° (cell), 

hK{g)= / n K (x)g(x)dx , g € L°° (cell) , (70) 

J cell 

then \tik\k converges weakly in L°° (cell)* . 

Proof. From Banach-Alaoglu theorem 120] one knows that the closed 
balls in L°° (cell)* are compact in the weak topology. Then, because 
\\nR\\ = No, it follows that the sequence {uk}k nas a ^ l eas t one accu- 
mulation point. Due to the fact that 

\\n K+1 - h K \\ = \\n K +i - n K \\ LHceU) — > 0, (71) 
we can conclude that there is one and only one accumulation point. ■ 

5 Application: The Local Density Approxi- 
mation 

In the local density approximation, the effective potential becomes: 

V [n] = (n-n Q )*v + v xc (n) , (72) 

where v xc is a function of n instead of a functional. The value of v xc (n) is 
equal to the exchange-correlation energy per particle of the corresponding 
infinite, homogeneous system. Thus, there will be no volume dependence 
for the exchange-correlation potential. The Hartree potential however will 
depend, in general, on the volume. For finite range interactions, the Hartree 
potential does not depend on the volume when the volume becomes larger 
than the range of the interaction. This case has been considered in Ref. 
jll| when the thermodynamic limit of the Hartree model was analyzed. We 
impose the following conditions on the two-body interaction and exchange- 
correlation potential. 

(PI) The singularity of the two-body potential is at least L 2 integrable. 
(P2) v (x) ~ \x\~ r as |x| — * oo with r > 2. 
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(P3) v xc : (0, oo) — > i? is diffcrcntiable and there exists p 1 such that 
t~ 1/p v xc (t) and t 1 " 1 /?— ^£ (f) are uniformly bounded over [0, oo). 

The conditions we imposed on the exchange-correlation include, for 
example, the case of an homogeneous electron gas. For this particular 
system, it was shown |2J that the low and high density behavior of the 
exchange-correlation potential is dominated by the exchange part which 
is proportional to n 1 / 3 . Thus, we can choose p = 3 in (PS) to include 
this case. Let us mention that, because of this behavior, the difference 
v X c (ni) — v xc (ri2) decays much slower than n\ — n-2 in the low density 
limit. Thus, the condition (C2) fails for this particular potential unless 
we can prove that the density of particles is larger than a certain strictly 
positive value. 

5.1 Estimates on the Hartree potential 

Let us point out that it is the Hartree potential that forces on us to consider 
only short range interactions. If one compares the conditions (P1)-(P3) 
with the conditions from Ref. |1U| . one can see an improvement because 
now we allow the interaction to decay as \x\ r with r > 2 instead of 3. The 
neutrality condition will play an essential role here. Unfortunately, we are 
still unable to include the Coulomb interaction in our theory. The reason 
is that the Hartree potential increases too fast as the system approaches 
the thermodynamic limit for long range interactions. 

Theorem 10 Suppose the conditions (PI) and (P2) are satisfied. Then, 

for n, n' G Sp er and 8 e [0, 2ir) 3 the following are true: 

i) There exists 7^, independent of K or 9 G [0, 2tt) 3 , such that: 




(73) 



ii) There exists 



L , independent of K or 9 £ [0,27r) 3 , such that: 



\\(n — n) * v\ 



< L H \\n-n'\ 



L 1 (ceJ!) - 



(74) 



Hi) The Hartree potential satisfies the condition (CS). 



Proof. Let us divide the Hartree potential in two parts, 




(75) 
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where vol^ contains the unit cell plus the adjacent cells and vol^ = 
T\vol^l For n, n' € Sg, r , it follows 



V { h ] [n\-Vg> [n 1 ] 



L 1 (ceiJ) 



(76) 



s$ TV sup / dx\v{\x,y\)\\\n-n'\\ L1{c 

ydvolt 1 ) J cell 



ill) 



where N is the number of cells in vol^ . The estimate makes sense because, 
if the singularity of the interacting potential is L 2 integrable, then it is also 
L 1 integrable. Moreover, 



[n] 



ro, ,n ^ 27ViVoWsup / dx \v(yi,x)v(\x,y 2 \)\, 

L 2 (cell) U J^j; 



where the supremum goes over all y\ and 2/2 € vol^\ Again, the estimate 
makes sense because the singularity of the interaction is L 2 integrable. The 
last inequality combined with 



L°°(cell) 



from Ref. JT], leads to: 

V ( h ] N (-iA^ + a)- 1 / 



^k\' 2 H/ll 



L 2 (cell) 



V™ [n] 



L 2 (cell) 



L 2 (cell) 

-1 



(77) 



(78) 



L°°(cell) 



< 2k x J 2 MN . 



sup / dx\v(\x,y\)\ ||/|| i2(ceH) 

j/^vol^ •* cell 



For the second term we write 

4 2) [n] (x ) - y^ 2) [n'\ (x) 
x. y + R 



(79) 



= E / - 

- £/ K 

- J cell L V 



x,y + R 



(n (y) - n (y)) dy 

x,R ) (n(y) -n' (y))dy, 



where the sum goes over all the sites of the crystal less the origin and its 
first neighbors. It is this place where we used the neutrality condition. 
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Then: 



4 2) [n] K] (x) 



(2) U'l 



(80) 



< ^2 sup u ( 



- x,y£cell 
R 



x,y + R 



x, R 



in — n 



^ 2^ sup u ( R — x + y) — v [R — x 



- x.y£cell 
R 



\L 1 (cell) ' 



where the last sum goes over an infinite lattice. Denoting £ = R/R, R 
\R\, it follows from condition (P2) that, for large R, 



v ( R — x + yj — v yR — x 
t-{x-y)/R 



(81) 



= ct.R- r 
€ ct.R- r 



i-x/R 



C-(x-y)/R - Z-2/R 



and using 



\a + b\ — \a ~ b\ ^ 2\b\ we can conclude: 



v (l 



R — x + y) — v [R — x 



< ct.R 



-r-1 



(82) 



Thus the sum in Eq. I|80l) converges and we proved that there exists Ljj , 
independent of K, such that: 



4 2) [n]-V® [n'\ 



<LP\\n-n'\\ L1(ceU) 



L^(cell) 



«S v ce aL§ \\n - n'\\ L1{ 



which automatically leads to: 
] [n]-V^[n']\ 

L x {cell) 



cell) 



and 



(83) 

(84) 
(85) 



This ends the proof of point i) and »). For point in), we notice that the 
difference between Vjy [n] and vj^ +1 ^ [n] is given by: 



E 

R£dA K 



and the sum contains a number of terms proportional to [K + 1) and R 
is proportional to K. Then the affirmation follows from Eq. (ISi't .B 



/ 


[„( 


x,y + R 


)-.( 




)] 


/cell 
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5.2 Estimates on the exchange-correlation potential 

Condition (F3) automatically leads to the condition (CI) for the exchange- 
correlation potential. 



Proposition 11 For n e Sp er , 



v xc (n) + a) 1 



7a 



Proof. 



\v xc (n) 



I LP (cell) 



cell 



-1/ ^ P -, 

n(x) v xc (n(x)) n(x)dx 



i/p 



su V t-^v xc (t)\\n\\% 



sC Nl /P sup t-V'Vvc (t) 



and using Eq. (|77ll 



v xc (n) 



iA s + a) 



f 



L 2 (cell) 

< \\vxc{n)\\ LP{cM) (-|A ? + a) -1 / 



1/2-1/p 
cell 



hxc {n)\\ LP 



LP (cell) 



L"(cell) 

(-lAj + a)" 1 / 



L°°(cell) 



(86) 



(87) 



v ce« ll^c V'VllLP(ce/;) IIJ llL 2 (ce/0 

Combining the above result with the results from the previous section, it 
follows that the condition (CI) is satisfied for the local density approxima- 
tion of the effective potential. We can also prove directly that the map T 
is continuous on 



Proposition 12 Let p ^ 1. Then 

\\v XC - V XC (7l2)|| 



LP(cell) 



^ p sup 



<cW ll n l-"2|lii P (ceH) , 



and consequently: 

||2> [ni] -T K [n 2 ] Hi^r) 
< 2\j$l h \\m - 7i 2 || L i (r) 



3-3/p 



sup 



<c (*) 



(89) 



(90) 
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for iii 2 e S 



K 

per • 



Proof. For ti, t 2 S R+ 

v xc {h) - v xc (t 2 ) = 



i/p 



dv x 



-alt 1 /? 



< p sup 

teR+ 



< p sup 

for p ^ 1. We can continue: 

(m) - w xc («2)|| L p( ceH ) 



dtVp 
* 1_1/ V xc (t) 



1/p _ 1/p 



1*1 - *2 



1/p 



(91) 



cell 



p sup 



t 1_1/ Xc (*)| |«i (x) - n 2 (x)\ 1/p ^ dx^j 



i/p 



= p sup 



* 1 - 1/ X c (t) IK-n 2 ||^ 



Li{cell) 



An important consequence of the above result is that any L 9 norm with 
q < p of u xc (m) — w xc (n 2 ) is finite. In particular: 



\Vxc (ni) - w^c (n2)|| £ i 



L 1 (ce/i) 



(92) 



^ i-i/p 



<4 (*) 



m - n 2 



1/p 



We notice that the limits on the chemical potential were based only on 
the condition (CI). Then, following the steps of the proof of point ii), 
Theorem 3, and Proposition 4 we have successively: 



\\T K [ni]-T K [n 2 ]\\ L1(T) 



(93) 



< 2A7 $ (\\(ni - n 2 ) * v\\ Ll{T) + \\v xc (ni) - (^2)|| L i (r) ) 



Of course, this result is far from condition (C3) because, for realistic 
exchange-correlation potentials we must choose p > 1 above. The above 
result however, is the best estimate one can get on the map T if only L 1 
estimates on the density of particles are used. To complete our analysis we 
need L°° estimates on the density of particles. 
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Lemma 13 Let A and W be two self-adjoint operators on L 2 (cell) such 
that A" 1 exists, 

sup (x, y)\ < x < °°> (94) 

and AW A is bounded. Then: 

|VF(x,x)| O ce «x 2 pWM||. (95) 
Proof. From the above conditions it follows that 

f (y) — A -1 (y, x) e L 2 (cell) . (96) 

Thus: 

\W(x,x)\ = (f,(AWA)f)^\\AWA\\\\ff L2{cell) (97) 
< v c eiix 2 \\AWA\\ ■ 

Theorem 14 For n € 5^ r , 

n-min - ct.A < 1 1 Tif NIIl« ^ "max + d.X, (98) 

where: 

n mi „ = — inf Tr$ FD (-^A s -n m ) (99) 

Vcell ee[0,27r) 3 

1 



sup Tr^ FD - /i M ) . 

«ceH ee[0,27r)3 

and ci. is if independent. 

Proof. In the previous Lemma, we choose: 

VF q = 4> FC (H^ - Mn ) - $fd (-|A 3q - Mn ) (100) 

and 

2 



^q=(-2- A « q +«) • (101) 

In this case A~ x can be computed exactly^!] and 

l^- 1 (x,y)\ < — L= J2 e-M*-t-K\. (1Q2) 
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Thus, the conditions on Aq in the previous Lemma are satisfied. Moreover, 
because 



W« = | ( 
+ ±e^ ( 



e 2 "r, — e 2 



cosh 



2 A » q )~j cosh 
-l 



(itf^-M*)]" 1 (103) 



£ (#^ q) - Mn)] " - cosh [f (-±A ?q - M 



2 V ™ 



and 



and 



cosh 



are of trace class, it follows: 



(104) 
(105) 

(106) 



1 

* 2 



cosh 



e 2 



2 i^n ^ 



4 ?, 



cosh 



cosh 



(-| A 5 q -M n ) 



Using Proposition 5, we can continue: 



(107) 



2 1 (! - A 7a) 4 ^6[eoToo) 



+7 F2 sup (z + a) e 2 (* m„) 

a;e[0,oo) 



sup (x + a) cosh § (a; — M„) 
0*) [, 



where Fi (z) = e 2 z and F 2 (z) = cosh ^ (z — jj, n ) 
V™ [n] 

L^ceH) 

< 2iY L ff + « c 1 e - 1/p ||^ c (n)|| iP(ceH) 



L!(ceH) 



. Moreover, 



^ 2iV L ff + iV 1/p ^ 1/p sup i" 1 ^ (t) . 
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Because fi n € [Mm>A*M]> we can conclude that 

PqW q A q || < pqWqAqH, ^ ct.X, 

where ct. is independent of q or K. Consequently, || Wq|lx,oo( ce ;n 
where ct. is again independent of q or K. Then the affirmation 
from 

T K [n\- (l + e^-^-^y 1 (x,x) 



(109) 



L°°(T) 

< ci.A 



L°°(ceH) 



and 



„ 'FD (-|A-/i n ) (x,f)|| ioo(r) 

- \\ t k [n] - <$>fd (-iA-/x„) (x,f)|| LOO(r) 

+ \\T K [n] - $ FD (-±A - (z, 
by observing that 

"min < \\$FD (-5A - fl n ) (x,x)\ 



(110) 



l£°°(T) ' 



I £00(3^ ^ ^max- 



(in) 



Let us now return to the last condition (C2) which remains to verified. 
We define the set B as a the strip in L°° (cell): 

B = {ne L°° (cell) , n min - e < n < n max + e, a.e.} , (112) 

where e is a positive constant such that n m i n > e. Observing that v' xc (t) 
is bounded over I — [n m ; n — e, n max + e], it follows: 

\\vxc (ni)-v xc (n 2 )\\ L1(cell) 
n2{S) dv xc (t) 

ni(x) 



cell 



dt 



dx 



(113) 



/ sup \v' xc (t)\ |ni (x) — 112 (x)\ dx 

J cell tei 



< su pKc 0)1 IK - n i\ 



te/ 



L!(ce«) 



where we omitted a set of zero measure where n (x) can have values which 
are not in the interval I. For A smaller than a certain value, independent 
of K, it follows from Theorem 14 that Tr- [S^ r ] C -Bk which completes 
our analysis of the local density approximation. 
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